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Based on 106 xlO 6 (3686) events collected with the BESIII detector at the BEPCII facility, a 
partial wave analysis of 0(3686) — ► ppn° is performed. The branching fraction of this channel 
has been determined to be 5(^(3686) -> pprv ) = (1.65 ± 0.03 ± 0.15) x 10" 4 . In this decay, 
7 N* intermediate resonances are observed. Among these ./V* resonances, two new resonances 
are significant, one 1/2+ resonance with a mass of 2300tf o t 1 09 MeV/c 2 and width of 340;^ t^ 
MeV/c 2 , and one 5/2" resonance with a mass of 2570t 1 9 l?o MeV/c 2 and width of 250t2 4 t2? 
MeV/c 2 . For the remaining 5 iV* intermediate resonances, the analysis yields mass and width 
values which are consistent with those from established resonances. 

PACS numbers: 14.20. Gk, 14.40.Lb, 11.80.Et 



Although symmetric non-relativistic three-quark 
models of baryons are quite successful in interpreting 
low-lying excited baryon resonances, they tend to predict 
far more excited states than are found experimentally 
("missing resonance problem") [l|, From the 

theoretical point of view, this could be due to a wrong 
choice of the degrees of freedom, and models considering 
di-quarks have been proposed Q. Experimentally, the 
situation is very complicated due to the large number 
of broad and overlapping states that are observed. 
Moreover, in traditional studies using tagged photons 
or pion beams P4llT|. both isospin 1/2 and isospin 
3/2 resonances arc excited, further complicating the 
analysis. 

An alternative method to investigate nucleon 
resonances employs decays of charmonium states such as 



J ftp and -0(3686). By selecting specific decay channels, 
such as -0(3686) — > ppir , N* intermediate resonances 
coupling to pir° or pir° can be studied. Here, A 
resonances are suppressed due to isospin conservation. 
As a consequence, the reduced number of states greatly 
facilitates the analysis (l2| . 

TV* production in J / yb — > ppr] was studied using partial 
wave analysis at BES [13[, and two N* resonances were 
observed. In a recent analysis of J/tp — > pnn~+c.c. [3], a 
new TV* resonance around 2000 MeV/c 2 named A r (2065) 
was observed. This A r (2065) was also observed in the 
decay of J/ip ppn pi]. The production of N(2065) 
in J ftp decays occurs close to the edge of the phase space. 
Thus, a similar search for this resonance in the 0(3686) 
decays should provide further insight. 

In the work of the CLEO Collaboration [3, 
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0(3686) -> ppTT was studied using 24.5 x 1O 6 0(3686) 
events. With the invariant mass spectra of pw 
and pp, two N* resonances (JV(1440), V(2300)) and 
two pp resonances(i?i(2100), i?2(2900)) were investigated 
without taking into account possible interferences 
between the resonances. The inclusion of R(2100) is 
suggested by a threshold enhancement in the pp mass 
spectrum. The concentration of events below 1800 
MeV/c 2 in the pir° mass spectrum is considered as the 
contribution of 7V(1440) alone. 

In this Letter, we report a study of N* resonances from 
■0(3686) —} ppir based on a data sample corresponding to 
106 million -0(3686) decays collected with the upgraded 
Beijing Spectrometer (BESIII), located at the Beijing 
Electron-Positron Collider (BEPCII) [H. 

The BESHI detector is composed of a helium-gas 
based drift chamber (MDC), a time-of-flight (TOF) 
system, a Csl (Tl) electro- magnetic calorimeter (EMC), 
a super-conducting solenoid magnet and a resistive 
plate chambers based muon chamber. More detailed 
information about the detector can be found in Ref . fl9| . 

The final state in this decay is characterized by two 
charged tracks and two photons. Two charged tracks 
with opposite charge are required. Each track is required 
to have its point of closest approach to the beam axis 
within ±20 cm of the interaction point in the beam 
direction and within 2 cm of the beam axis in the plane 
perpendicular to the beam. The polar angle of the track 
is required to be within the region of | cos(0)| <0.8. 

The TOF and the specific energy loss dE/dx of a 
particle measured in the MDC are combined to calculate 
particle identification (PID) probabilities for pion, kaon 
and proton hypotheses. For each track, the particle 
type yielding the largest probability is assigned. In this 
analysis, one charged track is required to be identified as 
a proton and the other one as an anti-proton. 

Photon candidates are selected by requiring a 
minimum energy deposition of 25 McV in the barrel 
EMC or 50 McV in the endcap EMC. To reject photons 
due to charged particle radiation production, the angle 
between the photon candidate and the proton is required 
to be greater than 10°. A more stringent cut of 30° 
between the photon candidate and anti-proton is applied 
to exclude the large number of photons from anti-proton 
annihilation. 

For events with one proton, one anti-proton and at 
least two photons, a kinematic fit (4C) with the sum of 
four-momenta of all particles constrained to the energy 
and three momentum-components of the initial e + e _ 
system is applied. A further kinematic fit (5C) with one 
more constraint of ir° mass for the two photons is applied 
to provide more accurate momentum information of the 
final states. When more than two photons arc found 
in a candidate event, all possible p|>77 combinations 
are considered and the one yielding the smallest Xbc * s 
retained for further analysis. 

The events passing the above selection criteria are 
shown in plots (a) and (b) of Fig. [TJ displayed as the 



Dalitz plot of 0(3686) — > ppir and the invariant mass of 
pp. The pp mass spectrum shows a clear J/0 signal. Due 
to the detector resolution, the observed width of J/0 is 
far larger than its natural width. This width difference 
causes a problem in the inclusion of J/0 in partial wave 
analysis. Thus, a cut of \M pP - Mj^\ > 40 McV/c 2 is 
applied to exclude events with pp arising from J/ decay. 
4988 events survive the event selection criteria. The mass 
spectra of pir° and pir° for the surviving events are shown 
in Fig. 12(c) and (d). 

For this analysis, two background sources are studied. 
The first one arises from 0(3686) decays and has been 
studied with two methods. In the first method, a 
sample of 10 s MC-simulated 0(3686) events is used and 
forty events survive the event selection, mainly due to 
misidentified or lost photons. In the second method, 
the background contribution is estimated using the ir° 
sideband events, defined by 30 MeV/c 2 < |M 77 -0.135| < 
45 MeV/c 2 . Only 26 events are found in the sideband 
area. The other background source arises from the 
continuum process e + e~ — > 7* — > ppir . This has 
been studied using 42 pb~ x of continuum data at y/s = 
3650 McV. After normalized to the integrated luminosity 
of 0(3686), 447 background events are found. In 
Fig. [2 the shaded histograms show the total background 
contributions from continuum process and ir° sideband, 
in which the continuum contribution accounts for about 
95%. 

In our present investigation, with larger statistics than 
at CLEOc, more than one N* state below 1700 MeV/c 2 is 
seen in the pir° and pir° mass spectra, and the threshold 
enhancement in the pp mass spectrum is also visible. 
To better understand the components of this decay, a 
partial wave analysis taking into account the possible 
interferences is pursued. 

The decay of 0(3686) —> ppir is supposed to be 
dominated by two-body decays involving N* , N* states, 
which can be described by 0(3686) — > pN*(pN*), 
N*(N*) pir°(pir°). In addition, a process of the 
type 0(3686) — > Rtt° is considered, where R represents a 
hypothetical pp resonance. The data are fitted applying 
an unbinned maximum likelihood fit. The amplitudes 
(Ai) for all possible partial waves are constructed using 
the relativistic covariant tensor amplitude formalism (]~5l — 
With these amplitudes, the total transition 
probability for each event is obtained from a linear 
combination of these partial wave amplitudes as u> = 
\T,iCiAi\ 2 . Finally the likelihood function ln(L) is 
constructed as V*™ n ln( f jfeA e wl ), where n is the total 
number of events, £ is the four-momentum of p, p and 
7T°, w(£) the probability density for a single event to 
populate the phase space at £, and e(£) is the detection 
efficiency to detect one event with £. The free parameters 
Ci are determined by maximizing the likelihood function 
ln(L). For each N* state, the amplitude is parameterized 
with a Breit-Wigner function, in which the mass and 
width of the resonance are variables, as described in 
The background contributions from ir° sideband 
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FIG. 1: (a) The Dalitz plot of 0(3686) -¥ ppTv", the invariant 
mass spectra of (b) pp, (c) pit , and (d) piv . The dashed 
lines in (b) show the cut at the J/0 mass region. The 
crosses represent the experimental data, and the shaded 
histograms show the background from continuum process and 
7T° sideband. The histograms in solid line show the sum of 
MC prediction and the background. 



best fitting result are taken as the optimized values. 
Table [TJlists the optimized values for the seven N* states. 
Here, the first errors are statistical and the second ones 
are systematic. In this table, the first five N* resonances 
are consistent with the values in the Particle Data Book 
[20j . while the last two states can not be identified with 
7V(2100) or 7V(2200). However, the significance of these 
two states are 15<7 and 11.7 a. As a consequence, we 
label these two states as 7V(2300) and iV(2570), with J p 
assignment of l/2 + and 5/2 - , respectively. 

Using these eight significant resonances, the fit result 
agrees well with the data, as shown in Fig. [T] The \ 2 over 
the number of degree of freedom is 1.12. The contribution 
of each intermediate resonance including interference 
effects with other resonances are extracted and shown 
in Fig. [H Plot (a) shows the contributions of iV(1440), 
7V(1520), 7V(1535) and 7V(1650) in which we can see clear 
peaks and also tails at the high mass region from the 
interference effects. Plot (b) shows the contributions 
of 7V(940), JV(1720), 7V(2300) an d 7V(2570). For 
-/V(2300) and A r (2570), their peak positions are below the 
Breit-Wigner mean values reported in Table I because 
of the presence of interference contributions, as well as 
phase space and centrifugal barrier factors. 



and continuum processes are removed by subtracting the 
log-likelihood (ln(L)) values, as the log-likelihood value of 
data is the sum of that of signal and background events. 
Possible interference between continuum processes and 
-0(3686) decays is not considered. 

All ./V* resonances up to 2200 MeV/c 2 with spin up 
to 5/2, listed in the summary tables of the Particle 
Data Book [2(|, are considered in this analysis, such 
as the well established states, iV(1440) and iV(1520), 
and not well measured states, 7V(2090) an d JV(2100). 
Phase space decay and two speculative N* resonances, 
iV(1885) and iV(2065) are also considered. According to 
the framework of soft n meson theory [2l[, the off-shell 
decay process is needed in this channel. Thus, iV(940) 
with a mass of 940 MeV/c 2 and zero width is included. 
The 7V(940) represents a virtual proton, which could emit 
a 7T°. The Fcynman diagram of this process can be found 
in [l5j |. In total, nineteen intermediate resonances are 
considered. The contribution from N* with spin 7/2 or 
larger is expected to be small due to the suppression by 
the centrifugal barrier for an orbital angular momentum 
L>2. 

In our analysis, the first step is to select the significant 
resonances among all these resonances. The significance 
of each resonance is determined from the difference of 
the likelihood values of fits with and without the given 
resonance, accounting for the change of the number 
of parameters. Resonances with significance greater 
than 5cr are taken as significant ones and include 
A^(940) and seven N* resonances. The remaining 
insignificant resonances are removed and only considered 
when estimating the systematic errors. The mass and 
width of N* states are varied, and the values with the 



TABLE I: The optimized mass, width and significance (Sig.) 
of the seven significant iV* resonances. AS represents the 
change of the log likelihood value. ANdof is the change of 
the number of free parameters in the fit. In the second and 
third columns, the first error is statistical and the second is 
systematic. The names of the last two resonances, iV(2100) 
and AT (2200), have been change to iV(2300) and JV(2570) 
according to the optimized masses. 



Resonance 
iV(1440) 
iV(1520) 
7V(1535) 
JV(1650) 
7V(1720) 
7V(2300) 
JV(2570) 



M(MeV/c 2 ) r(MeV/c 2 ) AS AN dof Sig. 

1 •3on+ il + 2i ^TTT ; TTT 

ijyu_ 21 _ 30 

1S10+2 + 11 



1510±f+i: 
1650 +5 + 11 

1DOU_ 5 _ 30 

17001^0+32 
2300t«lJ 09 

nc: 7n +19+34 



o 40 +4B+VU 
**U_40_156 
.,,.+20+0 
llo -15-40 
19n +20+0 
lzu -20-42 

150+ 21+14 
^OU_22_50 

/,c: n +109 + 149 

* JU _94 -44 

d4u -30-58 
9 c n +14+69 



72.5 
19.8 
49.4 
82.1 
55.6 
120.7 
78.9 



11.5a 



4 
4 
6 
4 



5.0o- 
9.3cr 
12.2cr 
9.6cr 
15.0cr 

11.7(7 



Various checks have been performed to test the 
reliability of this analysis. The first one is the spin 
parity check, in which the spin parity of each state 
of the optimized solution is changed to other possible 
values to test the other J p assignments. For 7V(2300) 
and ,/V(2570), 1/2+ and 5/2~, respectively, are the best 
J p values. The significance becomes worse by at least 
3.5cr using other J p assignments. The second one is the 
Input-Output check. A MC sample was generated with 
given components. After the fitting procedure described 
above, the significant states and their properties (mass, 
width, branching fraction, and the effect of interference 
terms) are compared with the input values. The output 
values agree with the input within ±lcr, corroborating 
that the analysis procedure is reliable. 





FIG. 2: The contribution of each intermediate resonance in 
the pn° mass spectra. The interferences with other resonances 
are included. 



determined as follows, 



E(V>(3686) -> ppir ) 



N - N bkg 



e x N, 



</>(3686) 



x B(ir° — > 77) 



= (1.65 ±0.03 ±0.15) x 10~ 4 

Here, N represents the number of observed events, 
Nf,kg stands for the number of estimated background 
events, and e is the efficiency derived from MC events 
generated according to the model derived from the PWA 
analysis. This result is in agreement with the value of 
(1.33 ±0.17) x 10" 4 in the Particle Data Book [H. The 
products of the production and decay branching fractions 
for each N* intermediate resonance are also determined, 
as shown in Table HU 



TABLE II: Summary of measurements of the number of 
events, the MC efficiency (e), and the branching fraction 
(B.F.) of each intermediate resonance and the whole channel. 
Here, for the number of events and the branching fraction, 
the first error is statistical and the second is systematic. 



Resonance 



^(940) 
AT(1440) 
AT(1520) 
AT(1535) 
iV(1650) 
AT(1720) 
iV(2300) 
iV(2570) 
Total 



1870 
1060 
190 
673 
1080 
510 
948 
795 



N 

+90+487 



<%) 



90-327 
+90+459 
90-227 
+ 14+64 
14-48 
+45+263 
45-256 
+77+382 
77-467 
+27+50 

27-197 
+68+394 
68-213 
+45+127 
45-83 

4515±93 



27.5 ±0.4 
27.9 ±0.4 
28.0 ±0.4 

25.8 ±0.4 
27.2 ±0.4 

26.9 ±0.4 

34.2 ±0.4 

35.3 ±0.4 
25.8 ± 0.4 



B.F.(xlO" 5 ) 
f. 49 +U.20 + l.V8 
D-^-0.20-1.28 
o t: H +0.25+1.59 
O.OO_ ,25—0 .84 
n C/i+0.05+0.22 
u -°^-0.05-0.17 
,+0.28+0.99 
-0.28-0.97 
,+0.28+1.37 
-0.28-1.66 
3+0.10+0.24 
-0.10-0.71 
,+0.28 + 1.12 
-0.28-0.64 
,+0.08+0.40 
-0.08-0.30 

16.5 ±0.3 ±1.5 



2.47Z 
3.76j 
1.79j 
2.62j 
2.13j 



On the basis of the eight significant states, a scan 
for additional resonances has been performed with 
different spin parity mass and width combinations. No 
extra resonance has been found to be significant. For 
7V(1885), the obtained significance ranges from la to 
1.2a depending on the mass and width. The largest 
significance is obtained at a mass of 1930 MeV/c 2 and 
width of 150 MeV/c 2 . The significance for 7V(2065) 
varies between 3.2er and 4c, where the maximum is 
obtained at a mass of 2140 MeV/c 2 and width of 250 
MeV/c 2 . We consider neither resonance as significant 
and do not claim any evidence. Besides the known and 
speculative ./V* resonances, a 1 pp resonance candidate 
described by the Breit-Wigner function has been added, 
as suggested by the near threshold enhancement in the 
pp mass distribution. Varying the width from 50 MeV/c 2 
to 300 MeV/c 2 and mass from 1800 MeV/c 2 to 3000 
MeV/c 2 with the step size of 10 MeV/c 2 , the largest 
significance obtained is 4cr at a mass of 2000 MeV/c 2 and 
width of 50 MeV/c 2 , indicating that no pp resonance is 
required to explain the threshold enhancement. 



The branching fraction of ^(3686) 



-+ 



ppil" 



IS 



The systematic uncertainty sources are divided into 
two categories. The first includes the systematic errors 
from the number of ?/>(3686) events (4%), MDC tracking 
(4% for two charged tracks), particle identification (2% 
for both proton and anti-proton), photon detection 
efficiency (2%), and kinematic fit (7%). These 
uncertainties are applicable to all branching fraction 
measurements. The total systematic error from these 
common sources is 9.4%. The second source concerns the 
fitting procedure, which includes the uncertainties from 
additional possible resonances, the uncertainties using 
different Breit-Wigner paramctcrizations for partial 
wave amplitude, the uncertainties from background 
estimation, the uncertainties from the J/ip exclusion 
cut, as well as the differences in the Input- Output 
check. These sources are applied to the mass, width and 
branching fraction measurements of intermediate states. 

In summary, we studied the intermediate resonances, 
including their masses, widths and spin parities, in the 
decay ?/>(3686) — > ppir . Two new N* resonances are 
observed, in addition to five well known TV* resonances. 
The masses and widths as well as the spin parities 
of the two new N* states have been measured. The 
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branching fractions of -0(3686) — > ppn and the product 
branching fractions through each intermediate N* state 
are measured. No clear evidence for Y(1885) or Y(2065) 
has been found. The hypothetical pp resonance has a 
significance of less than 4<r, indicating that the threshold 
enhancement most likely is due to interference of N* 
intermediate resonances. 
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